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Abstract 

We investigate general properties of Kerr-Schild (KS) metrics in n > 4 
spacetime dimensions. First, we show that the Weyl tensor is of type II 
or more special if the null KS vector k is geodetic (or, equivalently, if 
T ab k a k b — 0). We subsequently specialize to vacuum KS solutions, which 
naturally split into two families of non-expanding and expanding metrics. 
After demonstrating that non-expanding solutions are equivalent to the 
known class of vacuum Kundt solutions of Weyl type N, we analyze ex- 
panding solutions in detail. We show that they can only be of the type II 
or D, and we characterize optical properties of the multiple Weyl aligned 
null direction (WAND) k. In general, k has caustics corresponding to 
curvature singularities. In addition, it is generically shearing. Neverthe- 
less, we arrive at a possible 'weak' n > 4 extension of the Goldberg-Sachs 
theorem, limited to the KS class, which matches previous conclusions for 
general type III/N solutions. In passing, properties of Myers-Perry black 
holes and black rings related to our results are also briefly discussed. 

1 Introduction 

Kerr-Schild spacetimes [1] have played an important role in four-dimensional 
general relativity. In particular, all vacuum KS solutions have been known 
for some time [1-3]. They are a subset of algebraically special spacetimes, 
and the corresponding KS null congruence is a geodetic, shearfree, re- 
peated principal null direction of the Weyl tensor. Notably, they include 
the Kerr metric, arguably one of the physically most important known 
exact solutions of Einstein's equations in vacuum. Other well-known vac- 
uum solutions which can be put in the KS form are, for example, pp- 
waves and type N Kundt waves. In addition, the KS class can also admit 
(aligned) electromagnetic or matter fields, so as to include, for instance, 
the Kerr-Newmann metric, the Vaidya solution and pp -waves coupled to 
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a null Maxwell field or to pure radiation. (We refer the reader to [4] for a 
comprehensive review and for a number of original references.) 

In recent years, gravity in more than four spacetime dimensions has 
become an active area of ongoing studies, mainly motivated by the in- 
creasing interest in string theory and extra-dimensional scenarios. In this 
context, the KS ansatz led to the remarkable discovery of rotating vacuum 
black holes in an arbitrary dimension n > 4 [5] . However, several features 
of higher-dimensional gravity have proven to be substantially different 
from the four- dimensional case, and more insight would be desirable. In 
this paper we explore systematically general properties of KS metrics in 
higher dimensions. The main focus will be on Einstein's theory in vacuum, 
bearing in mind, however, that some of our results may have applications 
also in different contexts. 

In addition to the intrinsic interest for possible specific KS solutions, 
an advantage of the KS ansatz is its mathematical tractability, which 
enables one to study it analytically in full generality. We thus present it 
also to illustrate a concrete application of the recently developed higher- 
dimensional classification of the Weyl tensor [6, 7] and Newman-Penrose 
formalism [8,9] (cf [10] for another very recent application). This will 
also enable us to present a partial extension (restricted to KS solutions) 
of the Goldberg-Sachs theorem [11,12] to higher dimensions, thus making 
contact with the previous results of [8] . The plan of the paper is as follows. 

In Sec. 2 we consider general KS spacetimes and study constraints on 
the compatible form of the energy-momentum tensor. In particular, wc 
show that the null KS vector fc is geodetic if and only if T a bk a k b — 
(similarly as in four dimensions [4]). 

In Sec. 3 we specialize to the case with a geodetic KS vector fc and 
show that in that case the spacetime is of Weyl type II or more special. 

In the rest of the paper we further restrict ourselves to vacuum solu- 
tions. These naturally split into a non-expanding and an expanding class. 
In Sec. 4 we study the non-expanding class and show that it is equivalent 
to type N Kundt vacuum solutions. As a consequence, we also observe 
that n > 4 pp-waves cannot always be cast in the KS form, as opposed to 
the four- dimensional case. 

The more complex expanding case is analyzed in Sec. 5. From a subset 
of the vacuum equations we obtain an optical constraint, i.e. a purely 
geometric condition on the KS null congruence fc in the flat background. 
This is then used to show that expanding vacuum KS spacetimes are 
necessarily of Weyl type II or D. By integration of the Ricci identities we 
subsequently determine the dependence of the Ricci rotation coefficients 
(in particular, of optical quantities) on an affine parameter r along the 
KS congruence. This enables us also to find the r-dependence of the KS 
function and thus to study general basic properties of KS geometries. In 
particular, we discuss a 'weak' extension of the Goldberg-Sachs theorem 
valid for vacuum KS spacetimes, which are generically shearing. We also 
demonstrate, in the 'generic' case, the presence of a curvature singularity 
at a caustic of fc. This is in particular true for all expanding non- twisting 
spacetimes, on which we comment as a special subset of KS solutions. 

We conclude in Sec. 6 with a summary and final remarks. In ap- 
pendix A frame components of the Riemann tensor are given for KS 
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spacetimcs with a geodetic k. Subsequent appendices B-D contain several 
proofs and technical details related to the results presented in the main 
text. 



1.1 Preliminaries 

We use the notation of [8] (see also [9]) and, in particular, we introduce a 



scalar products are £ a n a = 1 and m^ a m a = &%j (o, 6 = 0, 1, . . . , n — 1), 
and the spacetime metric reads g a t = 2£( a rib)+5ijma^m^ (sum over i, j). 
Directional derivatives along frame vectors will be denoted by D = £ a V a , 
A = n a X? a , Si = m' l ' a V a . The full set of the corresponding Ricci rotation 
coefficients is defined in [8,9]. In particular, we will often use the definition 
£a-,b = Lcdm^m^ (sum over c, d). In fact, it will be convenient to 
adapt the frame vector I to coincide with the KS null vector k. When 
k is geodetic, we will denote by r the corresponding affine parameter. 
Quantities that do not depend on r will be denoted by a subscript or 
superscript index (e.g., Ho, s?^, etc.). 

Let us also anticipate here that we will be using five types of indices 
with different ranges, namely 

a, b, c, . . . = 0, 1, . . . , n - 1, i, j, k, . . . = 2, . . . , n - 1, 

a, /?, . . . = 2p + 2, . . . ,m + 1, p, a, . . . = m + 2, . . . , n - 1, (1) 

fj,, v = 1, . . . ,p, with 0<2p<m<n — 2, 

where p and in are fixed integers defined later on. Einstein's summation 
convention is employed except for indices /i, v (for which summation will 
be indicated explicitly) , or when at least one of the repeated indices is in 
brackets (e.g., there will be no summation over i in DAij — —2s^Aij), 
unless (only in very few exceptional cases) stated otherwise. Note that 
for indices i,j, ... we will not distinguish between subscripts and super- 
scripts since there is no difference between covariant and contravariant 
components. 

When referring to equations presented in previous papers it will be 
sometimes convenient to denote them by the corresponding equation num- 
ber followed by the reference number, e.g. eq. (1 If, [9]). 

2 General case 

In this section we define KS spacetimes and we study their general prop- 
erties, without imposing Einstein's equations and without requiring the 
KS null congruence to be geodetic. 

2.1 Ansatz and geodetic condition 

We study an n-dimensional spacetime with a metric in the KS form 



null frame = 



n, m 




n— 1). The only non-zero 



gab = rjab - 2Hk a k b , 



(2) 
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with r/ab being the Minkowski metric diag( — 1, 1, 1), TL a scalar func- 
tion and k a a 1-form that is assumed to be null with respect to n ab , i.e. 
rf b k a k b = (n ab is defined as the inverse of n ab and rf b / g ac g bd Vcd)- 
From (2) it follows for the corresponding vector that k a = rf b kb = g ab kb, 
so that k a is null also with respect to g ab , and that the inverse metric has 
the form 

g ab = V ab + 2Hk a k b . (3) 
One also gets for the metric determinant 

g = n = -l. (4) 

Straightforwardly from the definition of the Christoffel symbols we 
obtain 

r"f, c = -(Hk a k b ), c - (Hk a k c ), b + r, ad (Hk b k c ), d + 2Hk a k d (Hk b k c ), d , (5) 
from which 

r a bc k b k c = o, r a bc k a k b = o, r a ab = o, (6) 

and consequently also 

ka.bk = k a ^bk , k - b k — k b k . (T) 

In particular, this implies that k a is geodetic in the flat geometry r\ ab iff 
it is geodetic in the full geometry g ab . We also have 

r a 6c fc a = (Hk b k c )- d k d , r a bc k b = -(Hk a k c ); d k d , (8) 

ka-.ck = k a ^c k =0, k a ; bc k -\-k a - b k ; C = k a , bc k -\~k a ^ b k , c = 0. (9) 
Since r a aft = 0, we can express the Ricci tensor as 

Rab = ^ c a b,c — r c di ,r d ac . (10) 

After employing the above formulas, Einstein's equations for the projected 
component R a bk a k b read 

R ab k a k b = 2Hg ab (k a . c k c )(k b . d k d ) = n T ab k a k b . (11) 

From now on, when using a null frame I, n, we will make the 

convenient choice 

£ = k, (12) 

adapted to the KS ansatz. 
Using eq. (11) and (see [8]) 

k a . b k b = L 10 ka + L i0 m a l) (13) 
(sum over i), we can formulate the following 

Proposition 1 The null vector k in the KS metric (2) of an arbitrary 
dimension is geodetic if and only if the energy-momentum tensor satisfies 
T ab k a k b = 0. 

Then, using frame indices, the condition of proposition 1 reads Too = 
(or, equivalently, i? o — 0), and it holds iff the null frame components of 
T ab do not include a term proportional to n a n b . This condition is of 
course satisfied in the case of vacuum spacetimes, also with a possible 
cosmological constant, or in the presence of matter fields aligned with the 
KS vector k, such as an aligned Maxwell field (defined by F ab k a ~ k b ) or 
aligned pure radiation (i.e., T ab ~ k a k b ). 
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2.2 Optics of the KS vector field 



It is interesting to discuss how the optical properties of k a in the two ge- 
ometries g a i, and n a b are related, i.e. to compare the Ricci rotation coeffi- 
cients [8,9] Lab and L ab defined with respect to g a b and r/ab, respectively. 
In order to do so one has to set up 'null' frames for both metrics. Let fc, 
n, be the frame for the full spacetime metric g a b as discussed above. 
Then, from eq. (2) it follows that rj a b = 2k/ a [n b \ + Tik b \] + Sijrria m b 
(sum over Thus, fc, n, m' 1 ' is now a convenient frame for the fiat 

metric r/ ab , provided one takes 

n a = n a + Tik a - (14) 

Now, from eqs. (8) and (13) one gets 

k a . b = Kb - (DH + 2HL w )k a k b - 2HL io m^k h) (15) 

(sum over i). Expanding k a]b on the frame fc, n, and k a ,b on the frame 
fc, n, m^' in terms of L ab and L a b, respectively, and using (14) one finds 
that L a b = L a b except for Lu — Lu — DTi — TiLio and Lu — Lu — TiLm. 
In particular, one has Lio = Ljo (so that fc is geodetic with respect to g a b 
iff it is geodetic with respect to r/ a b, as already observed in section 2), and 
the matrix 

La = k a]b m (i)a m (j)b = k a , b m (i)a m (j)b (16) 

is the same with respect to both metrics. The optical scalars shear, twist 
and expansion [8, 9] are thus unchanged 1 (note that when fc is geodetic 
this statement is independent of the particular choice (14) since Lij is 
then invariant under null rotations with fc fixed [9]). 

3 The case of a geodetic KS vector field 

In the rest of the paper we will assume that the null KS vector fc is 
geodetic and affinely parametrized (i.e., Lio = = Lio), which seems to 
be the simplest and physically most interesting case (recall proposition 1). 
The following definitions [8,9] will thus be useful: 

Sij = = (Tij + 95ij, Aij = 

8 = -^r^Sii, a 1 = OijOij, uj 2 = AijAij. (17) 

We shall refer to Sij , Oij and Aij as the expansion, shear and twist matri- 
ces, respectively, and to 9, a and ui as the corresponding scalars. 
Under the geodetic assumption, eq. (15) reduces to 

ka;b = ka.b - (DH)k a k b , (18) 

which will be employed in the following. Using k a b = r\ ac k Ctb (and, 
from now on, denoting ( ),d?? cd by ( )' c , e.g. k a ' b = ri bc k a , c etc.) one 
Gctsily finds ka hk ' — ka -fok 5 k^^bk — - ka'b k ' - ka 7 k c — k a ' k ■ c , 
k a ,bc k b — —k a ,b k b , c , k a ,bc fc' ' k c = and similar identities, also to be used 

1 Cf [13] in the case n = 4. 
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in following calculations. In addition, certain expressions will take a more 
compact form if we write them in terms of the Ricci rotation coefficients 
and of directional derivatives as 

k jd — La — Sii] ka. b k — Lij Lij , ka. b k — LijLji^ 

H, a k a = DH, H, ab k a k b = D 2 H. (19) 
From eq. (5) we thus get 

r e ab , e = (nk a k b ); e - (Hk a k e ) M - (nk b k e ), ae 

+ 2 [UD 2 U + {DH + HLn)DTC\ k a k b , (20) 
r e fb T f ae = 2 [(DH) 2 - 2W 2 w 2 ] k a k b , (21) 

The Ricci tensor then reads 2 

R ab = {Hk a k b ) t e e - (Hk a k e ) >be - (Hk b k e ) iae 

+ 2H[D 2 H + LuDH + 2Hu 2 ]k a k b . (22) 

Consequently, k is an eigenvector of the Ricci tensor, i.e. 

R ab k b = — \D 2 H + (n - 2)6DH + 2H<J\k a . (23) 

For certain applications it may be useful to observe that the mixed 
Ricci components, i.e. R\ = (Hk a k b )% - (Hk a k s ) ibs - (Hk b k s ) t a s , are 
linear in H [14, 15]. The Ricci scalar is thus also linear and reads 

R = -2 [D 2 H + 2(n - 2)8DH + H{n - 2)(n - 3)<9 2 + H(cu 2 - a 2 )] . (24) 

So far we have worked in Minkowski coordinates adapted to the flat 
background metric rjij . Note, however, that expressions for scalars, such 
as frame components of tensors, do not depend on that choice. 

3.1 Frame components of the Ricci tensor 

From eq. (22), one finds the non- vanishing frame components of the Ricci 
tensor 

Roi = —\D 2 H + (n - 2)6DH + 2Hu) 2 \, (25) 
Rio = 2HL lk L jk - 2[DH + (n - 2)6H]S lJ , (26) 

Jin = 8i{8iH) + (Na - 2HLu)DH + (4ii 3 - 2L n - Mji)SjH - L U AH 
+ 27Y^25jL[ij] +4LijL[ij] + LuLn — LuLu 

+ 2L m Mu - 2A ij N ij - 2Hlo 2 ^, (27) 
Ru = -Si(DH) + 2L {ll] DH + 2L tJ SjH - L jd 5iH 

+ 2H (SjAij + AijMkk — AkjMkj — LjjLu + SLijL^j] + LjiL^j^ . (28) 
Note that Roo = = Roi identically. 



2 This formula for R ab is equivalent to eq. (32.10) given in [4] for n = 4, up to rewriting 
partial derivatives as covariant ones. 
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3.2 Algebraic type of the Weyl tensor 

The full set of the frame components of the Riemann tensor is given in 
appendix A for any KS geometry with a geodetic, affinely parametrized 
KS vector. Thanks to RoiOj = Roioi = Roijk = and to Roo = Roi = 0, 
for the Weyl tensor we find immediately 

Comj = 0, Cokh = 0, Coijk = 0. (29) 

The Weyl tensor components with boost order 2 and 1 are thus identically 
zero, which enables us to conclude 

Proposition 2 Kerr-Schild spacetimes (2) in arbitrary dimension with a 
geodetic KS vector k are of Weyl type II (or more special), with h being 
a WAND of order of alignment > 1. 

Note that the KS null vector fc must indeed be geodetic for a wide 
class of matter fields, in particular in vacuum (cf proposition 1), so that 
proposition 2 applies in those cases. 

Another interesting result follows from the observation [16] that space- 
times (not necessarily of the KS class) which are either static or stationary 
with 'expansion' and 'reflection symmetry' can be only of the Weyl types 
G, Ii, D or O (see [16] for details and precise definitions). Taking the 
'intersection' of this family with the set of KS spacetimes considered in 
proposition 2, we can conclude that in arbitrary dimension n > 4 

Proposition 3 All static spacetimes of the KS class with a geodetic k 
are of Weyl type D or O. All stationary spacetimes of the KS class with 
'reflection symmetry' and with a geodetic, expanding k are of Weyl type 
D or O. In both cases, if the Weyl tensor is non-zero, k is a multiple 
WAND. 

As an immediate consequence, the higher-dimensional rotating black 
holes of Myers and Perry (indeed obtained using the KS ansatz [5]), are 
necessarily of Weyl type D in any dimension, as we anticipated in [16] 
(this had previously been demonstrated by an explicit calculation of the 
Weyl tensor in [17]). This also applies to uniform black strings, since 
adding flat dimensions to a KS metric clearly preserves the KS structure. 
On the other hand, proposition 2 implies that five-dimensional vacuum 
black rings [18] do not admit a KS representation, since they are of Weyl 
type Ii [19]. Let us also emphasize that proposition 3 is not restricted 
to vacuum KS solutions. For instance, it can also be used to conclude 
that static black holes with electric charge (and, possibly, a cosmological 
constant) [20] are also of Weyl type D. 

3.3 Vacuum equations 

Vacuum solutions must satisfy R a b = 0. We note from eqs. (25)-(28) that 
the Ricci tensor components Roi and Rij are simple and do not involve 
Ricci rotation coefficients other than Lij, which characterize the optical 
properties of k. It is thus natural to start from the corresponding vacuum 
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equations. Imposing Rij — gives 3 

(DlnH)Sij = L lk L ]k - (n - 2)0 Si j. (30) 

Contracting eq. (30) with S lJ we obtain 

(n - 2)0{D\nH) = L lk L lk - (n - 2) 2 6 2 

= a 2 + cu 2 - (n - 2)(n - 3)6 2 , (31) 

while its tracelree part (i.e., Rij — ^z^Sij — 0) is 4 

(DlnW)ffy = (a 2 l3 - ^- 2 a 2 S l3 ) - (A 2 a + ^u,%) 

+ 2a k(i A j)k - (n - 4)6a l3 . (32) 

Next, the equation R i = requires 

D 2 H + (n - 2)6DH + 2Hcu 2 = 0. (33) 

Using eq. (31), this can also be rewritten as D 2 H = [— a 2 — 3lu 2 + (n — 
2)(n - 3)6 2 ]H. 

Note that eq. (31) involves the lunction TL in a non-trivial way only 
when 9/0. It will thus be convenient to study non-expanding and 
expanding solutions separately in the following sections. The remaining 
vacuum equations seem to be of little help in a general study, and there 
is no need to write them down at this stage. 



4 Non-expanding vacuum solutions 

4.1 Einstein's equations 

The case = turns out to be somewhat special, since from eq. (31) we 
immediately get a = = ui, and thus 

Uj = 0, (34) 

and from this and eq. (33) 

D 2 U = 0, (35) 
i.e., Ti = Tio + Gar (recall that in our notation r denotes an affine pa- 
rameter along k and DTio = = DQo). The vacuum equations Rij = 
and Roi = are thus identically satisfied. Note, in particular, that the 
sole non-expanding condition implies that we are restricted to a subset of 
the higher- dimensional Kundt class of non-expanding, non-shearing and 
non-twisting vacuum solutions. 

Setting Lij = in the Ricci components (27) and (28), the remaining 
vacuum equations read 

5i(6iH) + Nnim + (4L XJ - 2L d i - M 3l )& 3 H 

+ 2H(2S l L [u] + 4L lt L [u] + L a La + 2L wl M ii ) = 0, (36) 

5i(DH) - 2L [a] DH = 0. (37) 

3 Hereafter, for brevity we shall write H~ 1 D'H = DlnH, where it is understood that InH 
should be replaced by In \H\ whenever H. < 0. 

4 For n = 4 this reduces to (DlnH)aij = 2a k ^Aj^ k . 
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4.2 Equivalence with Kundt solutions of type N 

4.2.1 Non-expanding KS solutions are of type N 

We have seen above that vacuum non-expanding KS solutions are Kundt 
spacetimes. In general, the higher-dimensional Kundt class is known to be 
of Weyl type II or more special, provided Roo = = Rot (thus in particular 
in vacuum) [9]. However, here we show that, in arbitrary dimension, 
vacuum KS spacetimes of the Kundt class are restricted to type N. In 
fact, it has already been shown in section 3 that the components of the 
Weyl tensor with boost weight 2 and 1 vanish. Using eqs. (34), (35) and 
(37) together with eqs. (A2)-(A6) we also find that all Weyl components 
with boost weight and —1 are zero (since in vacuum R a bcd = C a bcd)- 
The Weyl tensor is thus of type N, q.e.d.. 

4.2.2 Kundt solutions of type N are KS 

We have thus demonstrated that vacuum solutions of the KS class with a 
non-expanding KS vector are a subset of Kundt spacetimes of Weyl type 
N. We now show that the converse is also true, namely that vacuum Kundt 
solutions of Weyl type N are of the KS form, so that the two families of 
solutions coincide. 

Vacuum Kundt solutions of Weyl type N belong to the family of space- 
times with vanishing scalar invariants (VSI) [21,22]. We can thus begin 
with the higher-dimensional VSI metric [22, 23] 

ds 2 = 2dtt \dr + H(r,u,x k )du + Wi{r,u,x k )dx^ + 6 ij dx i dx j , (38) 

with i,j,k = 2,...,n — l, k a dx a = du being the multiple WAND. 

Similarly as for n — 4, higher-dimensional vacuum Kundt spacetimes 
of Weyl type N consist of two invariantly defined subfamilies [22, 23]: 
Kundt waves (with Lu = Lu ^ 0) and pp -waves (for which Lu = La = 
= Lu). It is convenient to discuss these subfamilies separately. 

The metric functions H and Wi for higher-dimensional type N Kundt 
waves are [22] 

W 2 = -2-^, 

x 2 

W s = x q B qs (u) + C s (u), (39) 
H = ^ + H°(u,x 1 , 

where s, q = 3, . . . , n — 1, B qs — —B sq and H must obey a field equation 
given in [22]. As shown in [22], metric (38), (39) is flat for 

1 n—1 

H°(u, x l ) = H { L = ^yZW? + x 2 F Q (u) + x 2 x l F,{u), (40) 

where Fo(u) and Fi(u) are arbitrary functions of u. Therefore metric (38), 
(39) is in the KS form for an arbitrary choice of H°(u, x l ), since it can 
always be rewritten as ds 2 = ds 2 lat + (H° — i/° at )dix 2 . 
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A similar argument can be used to show that also higher-dimensional 
type N pp -waves belong to the KS class - the corresponding metric func- 
tions Wi and i?n at entering (38) can be found in [22]. 5 

We can thus summarize the results of this section in the following 
(see theorem 32.6 of [4] for n = 4) 

Proposition 4 In arbitrary dimension n > 4, the Kerr-Schild vacuum 
spacetimes with a non-expanding KS congruence k coincide with the class 
of vacuum Kundt solutions of Weyl type N. 



4.2.3 KS does not include all -waves for n > 4 

A comment on some differences between n = 4 and n > 4 dimensional pp- 
waves is now in order. By a natural extension of the n = 4 terminology of 
[4], in any dimension pp -waves are defined as spacetimes (not necessarily 
of the KS form) admitting a covariantly constant null vector 1, i.e. £ a; b = 
0. It then follows directly from the definition of the Riemann tensor that 

Robed? = 0. (41) 

For four- dimensional vacuum spacetimes this is equivalent to the def- 
inition of the type N [4]. It is also known that, in addition to being of 
Weyl type N, all four-dimensional vacuum pp-waves can be cast in the 
KS form [4]. 

By contrast, in higher dimensions eq. (41) is only a necessary, but not 
sufficient condition for type N. It only says that the type is II (or more 
special) [19], and in fact for n > 4 one can easily construct explicit vacuum 
pp -waves of types III, II and D, as we now briefly discuss. 

Vacuum pp -wave metrics of Weyl type III can be obtained directly by 
specializing results of [22] to the vacuum case. One simple five-dimensional 
example is metric (38) with 

W 2 = 0, W 3 = h(u)x 2 x 4 , Wa = h{u)x 2 x 3 , (42) 



H = H = h(u) 



^((^y + ^VftV,* 3 ,* 4 ) 



(43) 



where h°(x 2 , x 3 , x 4 ) is 

In addition, one can also construct (ni + n2)-dimensional vacuum pp- 
waves of Weyl type II, e.g. by taking a direct product of a ni-dimensional 
vacuum pp -wave of Weyl type N or III with an Euclidean n2-dimensional 
Ricci-flat (but non-flat) metric (with both ni,ri2 > 4, and m > 5 if the 
ni-dimensional pp-wave is of type III). 6 Similarly, (m + n2)-dimensional 



5 There is a typo in eq. (96) of [22]: just drop the inequality m < n in the second sum. We 
thank Nicos Pelavas for correspondence on this point. 

°The proof of these statements is straightforward and we just sketch it. First, all the 
mentioned products arc really -waves since a covariantly constant vector field t defined 
in the ni-dimensional geometry can be trivially lifted to the (ni + ri2)-dimcnsional product 
geometry, in which it will still be covariantly constant. Such products are necessarily of 
Weyl type II or more special (with the lifted £ being a multiple WAND), as follows from the 
decomposable form of the Weyl tensor of product geometries [6, 16] . They cannot be of Weyl 
type III or N since they inherit non-zero curvature invariants (e.g., Kretschmann) from the ni- 
dimcnsional Euclidean space. They cannot be of type D since, starting from a frame adapted 



10 



vacuum pp -waves of Weyl type D arise if one takes a direct product of a 
m- dimensional fiat spacetime with an Euclidean n2-dimensional curved, 
Ricci-flat space (with m > 2, ri2 > 4). 

As a direct consequence of proposition 4, such n > 4 pp -waves of 
Weyl type II, D or III do not admit a KS form. They thus represent 
a 'counterexample' to both the n = 4 results, namely they are (higher- 
dimensional) vacuum pp -waves that are neither of type N nor of the KS 
form. 

Note, in addition, that pp -waves of Weyl type II and D cannot belong 
to the VSI class of spacetimes (which is compatible only with types III, 
N, O [21]), and therefore they necessarily possess some non- vanishing cur- 
vature invariants. In table 1 we summarize the aforementioned properties 
of various types of vacuum pp -waves in four (4D) and higher dimensions 
(HD). 7 





Weyl type 


KS 


VSI 


4D 


N 


V 


V 


HD 


N 


V 


V 


HD 


III 


X 


V 


HD 


11(D) 


X 


X 



Table 1: Properties of various types of vacuum pp -waves. 



5 Expanding vacuum solutions 

Non-expanding vacuum solutions of the KS class have been fully classified 
according to the discussion of the previous section. Let us now consider 
solutions with an expanding k, i.e. 6 =fc 0. 

5.1 Consequences of Einstein's equations 

For 9 / 0, from (31) one can write DhiTL as 

mnW=i^-(n-2)0. (44) 
Substituting into (30) one gets 

LikLjk = o\7i^ i i' (45) 



to the canonical form of the (type N or III) Weyl tensor of the n\ -dimensional spacetime, the 
components of negative boost weight turn out to be unchanged under null rotations, and thus 
cannot be set to zero. The only possible Weyl type is thus indeed II. 

7 In fact, in view of our previous comments a similar table applies to all Kundt solutions, 
in which case type III appears also in 4D [4], again VSI but not KS. 
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Interestingly, this equation is independent of the KS function Tt and it 
is thus a purely geometric condition on the KS null congruence k in the 
Minkowskian 'background' r\ a \, (recall the discussion of subsection 2.2). 
We will thus refer to it as the optical constraint. 8 Its consequences on the 
form of Lij will be studied in detail in appendix C and D, and they will 
be employed in the following sections. 

Using (17), eq. (45) can also be rewritten in terms of the shear and 
twist matrices as 

(n - 2)6 [(V, 3 - ^(7%) - (A 2 ,, + + 2a k(i A j)k ] 

= [a 2 + lo 2 - (n- 2)6 2 } cr i3 . (46) 

This traceless equation is in fact equivalent to (32), after using (44) (since 
the trace of eq. (45) is obviously an identity). 9 Contracting (46) with o~ij 
one finds the scalar constraint 

a 2 [a 2 +lo 2 - (n - 2)6 2 } - (n - 2)9 [a 2 ^^ - a i . j A 2 ji ] = 0. (47) 

So far we have discussed consequences of the vacuum equation Rij — 0. 
Next, one has to make sure that eq. (33) (i.e., Roi = 0) is now compatible 
with (44). In fact, taking the D-derivative of (44) and using the scalar 
Sachs equations of [9] and eq. (47), one exactly recovers eq. (33). This is 
thus automatically satisfied, provided the preceding equations hold. Of 
course, when looking for an explicit solution one should also solve the 
remaining Einstein equations, namely Rn — and Ru = 0. These are, 
however, too involved unless one makes some further assumptions (such 
as the presence of symmetries, etc.). In any case, they are not needed in 
the following discussion. 

5.2 Algebraic type 

From the general result of proposition 2, we already know that the Weyl 
tensor of any KS metric with a geodetic k is of type II or more special. 
Here we show that, in fact, the types III and N are not possible for vacuum 
expanding solutions. 

By reductio ad absurdum, let us thus assume that all components of 
the Riemann (Weyl) tensor with boost weigh zero, i.e. (A2)-(A4), vanish. 

8 This constraint is of course identically satisfied in the trivial case a = = lu (i.e., Lij = 
Sij = OSij) which includes, e.g., the Schwarzschild-Tanghcrlini solution [20]. Less trivial 
examples are provided by the KS congruence of static black strings, for which a ^ and lo = 0, 
or of rotating black strings and Myers-Perry black holes (which were indeed constructed in [5] 
using the KS anstatz), both having <r ^ ^ lo. (To explicitly verify these statements one may 
use the corresponding optical quantities calculated in any dimension in [8,16].) 

On the other hand, a simple example of a geodetic null congruence violating the optical 
constraint is provided by the vector field k = \/l + f 2 (4>)dt + f(4>)d p + d z in the flat geometry 
ds 2 = -dt 2 +dp 2 + p 2 d0 2 +dz 2 + S AB dy A dy B (with A, B = 1, . . . , n - 4). This geodetic null 
congruence is expanding and shearing for / ^ 0, and violates the optical constraint if / ^ ^ 0, 
in which case it is also twisting. 

9 In odd dimensions, it implies that one cannot have dij = if Aij ^ (this was already 
known in a more general context [9]). 
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In particular, we can multiply Ronj by Lij. Using the optical constraint 
(45), the condition Ron j Lij — gives 



SuDlnH 



(48) 



where we have dropped an overall factor Lj k Lj k = a 2 + uj 2 + [n — 
2)9 2 > 0. Using as a 'spatial' basis an eigenframe of Sij, we have 
Sij — diag(s( 2 ), S(3), S(4), . . .)• At least one eigenvalue must be non-zero, 
say S( 2 ) 7^ 0. Then from (48) with i = / = 2 we get 



S( 2 )D\nH = -2s ( 2)A 2 2, 



(49) 



which obviously can never be satisfied since A22 = (and since the case 
DH = is 'forbidden', cf appendix B). This contradiction completes the 
proof of 

Proposition 5 In arbitrary dimension n > 4, Kerr-Schild vacuum space- 
times with an expanding KS congruence k are of algebraic type II or D. 

For a similar result for n — 4 cf [4, 24] . 



5.3 Fixing the r-dependence 

As detailed in appendices C and D, the optical constraint (45) and the 
Sachs equations [9] imply that there exists an appropriate frame, satisfying 

i 

Mjo = 0, such that the matrix Lij takes the block diagonal form 

\ 



(i) 




(50) 



V 



The first p blocks are 2x2 and the last block C is a (n— 2— 2p) x (n— 2— 2p)- 
dimensional diagonal matrix. They are given by 



-(*«) 



S(2 M ) 



S(2„) 

r 

hT(^P 



j42/j,2/j+1 
S(2 M ) 



• ,P), 



A- 



2/j,2/j+1 



- 2 + K 2M) ) 2 ' 



£ = -diag(l, 
r 



1,0^^0), 

(m — 2p) (n — 2 — m) 



(51) 
(52) 



with 0<2p<m<n — 2. The integer m denotes the rank of L^, so that 
Lij is non-degenerate when m = n — 2. 
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Accordingly, the optical scalars are given by 

(n-2)e = 2± - T -^- + (53) 



, . , (a° 

M=l \ V (2m)'' 

+ (n-2-m)6> 2 . (55) 

Note, in particular, that expansion 8 indicates the presence of a caustic 
at r = 0, except when 2p = m (with m even). Note also that shear is 
generically non-zero (the special case a = will be discussed below). The 
twist is zero if and only if p — (or, equivalently, all oL\ vanish) . 

One can similarly fix the r-dependence of all Ricci rotation coefficients 
(at least with the additional 'gauge' condition iVio = 0). Since this is not 
needed in the following discussion, corresponding results are relegated to 
appendix D. 

Knowing the form of Lij enables us to solve the vacuum equation (44), 
where now L ik L ik = (n — 2)9r~ 1 . The r-dependence of Tt is thus given 
by 

This includes in particular the solution of (44) in the case when its rhs 
vanishes, i.e. DTL = 0, which happens for m = 1 (implying p = 0). This 
is, however, incompatible with the Bianchi identities [8] as explained in 
appendix B. Hence, in the following we shall restrict to 

2 < m < n - 2. (57) 

For example, the r-dependence in the special case of Myers-Perry so- 
lutions is obtained by setting m = n— 2 in (56) (note, however, that func- 
tion (56) is more readily compared with that of [5] using the parametriza- 
tion (6) of [25] for the direction cosines). Black holes with only one non- 
zero spin [5], in particular, correspond to the further choice p = 1, and 
static black holes to p = (see also subsection 5.7). 

In general, the asymptotic behavior of 7i for r — » oo is given by 

^=^k+0{r-^), (58) 

i.e., the function H behaves as a Newtonian potential in (m + 1) space 
dimensions with r as a radial coordinate. 



5.4 On the Goldberg-Sachs theorem 

Let us now use the above results to comment on a partial extension to 
higher dimensions (but restricted to KS solution) of the Goldberg-Sachs 
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theorem. In four dimensions, this is a well-known theorem stating that in 
a vacuum (non-flat) spacetime, a null congruence is geodetic and shearfree 
if and only if it is a multiple principal null direction of the Weyl tensor 
[4,11,12]. 

For n = 4, the matrix Lij associated with a generic null vector I [8] 
is 2 x 2 and can be written in terms of the standard Newman-Penrose 
notation [4] as [9] 10 

T --\( (p + p) + (° + °) -i(p-p)+i(c7-a)\ 
2\i{p-p)+i{a-a) (p + p)-(a + a) J' 

In vacuum, when I is a multiple principal null direction of the Weyl 
tensor the Goldberg-Sachs theorem implies a = 0, so that Lij reduces to 

13 2\i(p-p) p + p J 1 ' 

We have observed in previous sections that for vacuum KS spacetimes 
the KS vector I — k is indeed a (geodetic) multiple WAND for any n > 4. 
Accordingly, in the special case n = 4, eqs. (50)-(52) must be of the 
shearfree form (60). It is easy to verify that this is indeed the case, since 
for n = 4 we necessarily have m = 2 (cf (57)) and therefore Lij either 
coincides with one of the blocks £( M ) (for p = 1, i.e. p — p / 0) or 
is proportional to the two-dimensional identity matrix (for p — 0, i.e. 
p-p = 0). 

However, it has been pointed out [5, 8, 9, 16, 26] that the Goldberg- 
Sachs theorem cannot be extended to higher dimensions in the most di- 
rect formulation 'multiple WANDs are geodetic and shearfree in vacuum'. 
Indeed, we have already noticed that k is generically shearing. Neverthe- 
less, for any n > 4 our result (50)-(52) can be viewed as a 'weaker' higher- 
dimensional Goldberg-Sachs theorem, albeit restricted to KS spacetimes: 
in a suitable basis, the matrix Lij (associated with the geodetic, multiple 
WAND k) consists of 2 x 2 blocks that are 'shearfree', i.e. reflecting the 
four-dimensional shearfree condition (60) in various orthogonal 2-planes. 11 
The only 'exceptions' to this of course arise in odd spacetime dimensions 
(there will be an isolated one-dimensional block) and in the degenerate 
case det L — 0, in which one may add an arbitrary number of zeros along 
the diagonal (most simply, by just taking a direct product with flat dimen- 
sions). This 'generalized Goldberg-Sachs condition' can also be expressed 
in a basis-independent form simply as (dropping the matrix indices) 

[S,A]=0, A 2 = S 2 -FS, (61) 

where S = (L + L T )/2 and A — (L — L T )/2 are the symmetric and 
antisymmetric parts of L, and T can be fixed as in (CI) by taking the 
trace. The first equation in (61) is clearly equivalent to the condition 



10 Note that, exceptionally in this subsection only, the complex shear a does not coincide 
with the real shear scalar defined in section 3 and used throughout the paper (although they 
are simply related), but it is the usual Newman- Penrose spin coefficient. 

lx Let us recall that we have been focusing on the expanding case 9^0, since 9 = => 
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[L, L T ] — 0, i.e. L is a (real) normal matrix. The second equation in (61) 
is just a consequence of the first one when n — 4. 

Let us also emphasize that the above 'weak Goldberg-Sachs theorem' 
has, in fact, been proven in more generality (i.e., without the KS assump- 
tion) for higher- dimensional vacuum spacetimes of Weyl type III and N: 
the WANDs of such spacetimes must be geodetic, and the associated ma- 
trix Lij is indeed of the form (50) , with only one non-zero 2x2 block and 
zeros elsewhere [8]. 12 

For a possible complete extension of the Goldberg- Sachs theorem, one 
should study the remaining possibilities, i.e. n > 4 vacuum spacetimes of 
Weyl type II and D that do not belong to the KS family. In addition, one 
should also analyze the inverse implication, that is, whether in vacuum 
the existence of a geodetic null congruence with an associated matrix Lij 
of the form (50) implies that the algebraic type is II or more special (this 
is already known in the simple Kundt case Lij — [9]). Let us recall 
that the 'geodetic part' of the n > 4 Goldberg-Sachs theorem has been 
discussed in [16]. 

It is worth mentioning that a n > 4 extension of the Robinson theorem 
has been proven in even dimensions [27] assuming a generalization of the 
shearfree condition (see also [28-31]) different from ours. The relation to 
our work will be discussed elsewhere. 



5.5 Weyl tensor components of boost weight zero 

The r-dependence of the matrix Lij determines explicitly also the in- 
dependence of the Weyl components with boost weight (cf appendix 
A). In order to express these components compactly, it is convenient 
to introduce a (n — 2) x (n — 2) matrix <3>ij = Couj along with its 
symmetric and antisymmetric parts = &tij) and $^ = $\ij]- All 
boost order zero components of the Weyl tensor are then determined 
in terms of the (n — 2)(n — 3)/2 independent components of <E>^ and 
the (n — l)(n — 2) 2 (n — 3)/12 independent components of Cijki, since 
$g = -\Cikjk, Com = 2*£, Coioi = $ = [8, 16]. 

It follows from appendix A that, using the adapted frame of ap- 
pendix C, the matrix 3>ij inherits the block structure of Lij, with the 
only non-zero components and trace given by 

<E>2 M ,2 M = "I > 2 M +1,2 M +1 = — 27^2^2^+1 — S(2ii)DH, (62) 

$2 M ,2 M +i = *&, 2m +i = -D(WA 2 „, 2 „+i), (63) 

= --S a gDH, $ = D 2 U = -(n-2)9DH-2Hcu 2 . (64) 
r 

Up to index permutations, the non-zero Cijki components are 

C2 M ,2 M +1,2 M ,2 M +1 = 2H(3Al^ 2 ^+i ~ S (2 M )), (65) 
C2/i,2fi+l,2i/,2i/+l = 2C2 M , 2i/,2/i+l, 1v+\ = — 2C2^,2i/+l,2^+l,2i/ 

= 47M2 M ,2,x+i J 42 I /,2 I /+i, (66) 

12 More precisely, this has been proven for all n > 4 solutions of Weyl type N, for n > 4 
non-twisting solutions of type III, and for all n = 5 solutions of type III. For twisting type III 
solutions with n > 5 an extra assumption on the Weyl tensor was necessary [8]. 
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C2 m ,2i/,2/j,2i/ = C2 m ,2i/+1,2/j,2i/+1 = — 27iS(2 M )S(2„) , (67) 

C(a)(i)(a)(i) = —2Hs {i )r 1 , (68) 

where 

Let us observe that all the above components fall off at r — » oo as 
l/r m+1 or faster (in the non-twisting case, i.e. when all ^2^,2^+1 vanish, 
only l/r m+1 terms are present, cf also [10]). 

Since (as shown above) expanding KS spacetime can be only of Weyl 
type II or D in vacuum, Weyl components of boost weight fully determine 
the Kretschmann scalar, i.e. 

RabcdR abCd = 4(i?010l) 2 + RijklRijkl + 8RojliRoilj ~ 4i?01ij Roiij ■ (69) 

In vacuum R a bcd = C a bcd, and we can use the above compact notation to 
reexpress 

RabcdR abcd = 4$ 2 + djkiCijki + 8*g$g - 12 !<!>. '<!■,}. (70) 

This will be useful soon in the discussion of singularities. 

To conclude, let us observe that Weyl components of boost weight — 1 
and —2 contain also derivatives different from D, and we cannot study 
their r-dependence at this general level. 

5.5.1 Type D KS spacetimes: also the second multiple 
WAND is geodetic 

Type D KS spacetimes are a special subclass of general KS metrics. By 
definition of type D, only boost weight components of the Weyl tensor 
(as given above) are non-zero in an adapted frame, i.e. they admit a 
second multiple WAND not parallel to the KS (geodetic) vector fe. Here 
we show that this second WAND is also geodetic (as already known in the 
special case of Myers-Perry black holes [17], see also [16]). 

First, note that, from (63), one has <E>^ = <^> uj 2 = (cf eqs. (51) 
and (56)). In addition, since / (cf (64) and appendix B), for 
w = (i.e., p = 0) it is easy to see from (64) that / — $ for all 

values of i. Thus, for type D KS spacetimes proposition 6 of [16] implies 
that the second multiple WAND is indeed geodetic (since one can always 
align the frame vector n to it, and this is compatible with the assumed 

i 

condition Mjo = 0) . This is a general consequence of the Goldberg-Sachs 
theorem in four dimensions, but it is a non-trivial result for n > 4 (cf [16]). 
Therefore, it can be viewed as another partial extension of the Goldberg- 
Sachs theorem, again limited to KS spacetimes. 

5.6 Singularities 

The form (56) of the KS function suggests there may be singularities at 
r = 0, except in the cases 2p = m (m even) and 2p — m — 1 (m odd). In 
order to invariantly identify possible curvature singularities, let us analyze 
the behavior of the Kretschmann scalar (70). 

Eq. (70) clearly consists of a sum of squares, except for the last term, 
which is negative. For our purposes, it suffices to focus on the first and 
the last terms determined by (63) and (64) with (51) and (56). 
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5.6.1 'Generic' case (2p ^ m, 2p ^ m — 1) 

Excluding for now the special cases 2p — m and 2p = m — 1, from (56) 
one easily finds that for r ~ 

K~r- (m - 2p_1) , D-H~r" (m " 2p) , ~ r" (m " 2p+1) . (71) 

Inserting this into (63) and the second of eqs. (64), it is clear that 
the first term in (70) will dominate over the last term near r ~ 0, and 
therefore the Kretschmann scalar will diverge at r = 0, thus confirming the 
presence of a curvature singularity. Note that this 'generic' case includes, 
in particular, all non- twisting solutions, for which p = (and $^ = 0). 

As an explicit example, n > 5 Myers-Perry black holes with only one 
non-zero spin (m = n — 2, p = 1) fall in this subclass, and in appropriate 
coordinates [5] one has 2H. — —fir''~ n /(r 2 + a 2 cos 2 6) (a is the angular 
momentum parameter). See [5] for a detailed discussion, including other 
possible examples with more than one rotation parameter. 

5.6.2 Special cases 2p = m and 2p = m — 1 

Let us now comment on the cases 2p = m (m even) and 2p = m — 1 (m 
odd), for which H = rHo U™=i k 2 +(«(2 M )) 2 ] _1 andH = Ho Ilj l =i" 1)/ V + 
(a°2 M )) 2 ] _1 , respectively. With these assumptions TL and its D-derivatives 
are clearly non-singular at r — (since p denotes the number of non- 
vanishing terms fl?2 M ) m (56)). However, in general, fl?2 M ) may be functions 
of spacetime coordinates different from r. Given a specific KS solution, 
there may thus exist 'special points' where some of the a%ix) vanish. Let 
us say, e.g., that q (with q > 1) of such terms a?2 M ) vanish simultaneously 
there. Then, proceeding as above one finds that at those special points, 
when r — > 

H~r- 2g+1 , DH~r- 2q , D 2 H~r- 2q - 1 (2p = m), (72) 
H~r- 29 , Dn~r- 2q -\ D 2 n~r-' 2q -' 2 (2p = m - 1).(73) 

Again, the first term in (70) will dominate over the last one and the 
Kretschmann scalar will diverge at the 'special points' when r — > 0. 

For example, for both the n = 4 Kerr solution (m = 2, p = 1) and 
the n = 5 Myers-Perry metric with only one spin (m = 3, p = 1) the 
only non-zero «(2 M ) function is determined by (o^ 2 )) 2 ~ fl2 cos2 Hence, 
there is a curvature singularity at r = and 8 = ir/2 (this is the well- 
known ring-shaped singularity of the Kerr spacetime). By contrast, the 
n = 5 Myers-Perry solution with two non-zero spins (again, m = 3, p = 1) 
contains (d( 2 )) 2 = a 2 cos 2 8 + b 2 sin 2 (in the notation of, e.g., [26]), which 
never vanishes. Possible singularities of rotating black hole spacetimes are 
studied in more generality in [5] . 

5.7 Special subfamilies 

We now briefly comment on special subfamilies of expanding KS solu- 
tions, characterized by a non-twisting or a non-shearing KS vector k. 
In relation to the first possibility, let us note that general properties of 
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higher- dimensional non-twisting vacuum spacetimes (not necessarily KS) 
have been recently studied in [10], and that all such solutions are explic- 
itly known with the additional assumptions of non-zero expansion and 
vanishing shear ('Robinson- Trautman spacetimes') [32,33]. 

5.7.1 Non-twisting solutions 

The KS vector k is non-twisting if and only if Lij is a symmetric matrix, 
i.e. for p = 0. This clearly gives 

Lij = Sij = -diag(l, . . . , 1, 0, . . . , 0), (74) 



r 



(ra-2-m) 



so that (cf (53)-(55)) 



1 m „ 2 1 m(n — 2 — m) 

r n — 2 r 2 n — 2 

From eq. (56) , the r-dependence of H becomes simply 

n = ^j. (76) 

For Weyl components with boost weight we now have the non-zero 
components 

*»0 = <s = ;St("»-i)*«/s> ( 77 ) 

C (*)(«)W = - ^' ( 78 ) 

As mentioned above, all non-twisting KS solutions contain a curvature 
singularity at r — 0. 

Explicit n-dimensional examples of non-twisting KS spacetimes with 
m > 2 are static (uniform) black strings and black branes constructed as 
a direct product of a (m + 2)-dimensional Schwarzschild black hole with 
(n — m — 2)-dimcnsional flat space (which are of Weyl type D [16]). In 
suitable coordinates, these can be written as 

ds 2 = [-(2dr + KAu)Au + r 2 dtf m + 8 AB Ay A Ay B ~\ + ^-^du 2 , (79) 

where the metric in square brackets represents a flat n-dimensional space- 
time (with A,B — 1, . . . , n — m — 2, and dfi^ being the line element 
of a m-dimensional space of constant curvature K), k a dx a = du and 
2H = -fir- m+1 . 

The above non-twisting KS solutions are, in addition, non-shearing iff 
m — n — 2 (the Kundt case m = is ruled out here by the assumption 
6 7^ 0). In this case, they must belong to the family of higher-dimensional 
Robinson- Trautman spacetimes [32]. In fact, one can see from the results 
of [32] that the only Robinson- Trautman solutions that are also KS are 
given by static Schwarzschild- Tangherlini black holes, i.e. by metric (79) 
with m — n — 2 (in which case a cosmological constant can also be added 
[32] while still preserving the KS form). 13 



3 The n > 4 Robinson- Trautman family contains, in addition, generalized black holes with 
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5.7.2 Non-shearing solutions (n even) 

Non-shearing, non-twisting solutions belong to the already discussed KS- 
Robinson-Trautman class, so that we can now focus on non-shearing but 
twisting solutions. Recall [9] that a — and ui =fc can occur only for 
even n (and for 9 0). 

We see from (55) that a — requires m = n — 2 (i.e., Lij must be 
non-degenerate), so that m is also necessarily even. Then, we have to 
consider two possible situations. 

1. If 2p = n — 2, k is shearfree iff S( 2fJ ) = 9 for any fi — 1,2,..., (n — 2)/2 
(cf (55)), so that a® 2 ) = a(4) = • • • a(n-2) = a o, and 

ao 



V^a-TT^- (80) 



This agrees with the general behavior found in [9] (note also that 
AikAjk = -^r^Sij, as expected from the optical constraint (45) with 

(7 = 0). 

The function H reduces to 

n = — (si) 

This includes, for instance, the n = 4 Kerr solution. Explicit solu- 
tions for n > 4 seem to be presently unknown. 

2. If 2p < n — 2, we have the additional condition 9 — 1/r. This implies 
lo = ao — 0, which is the already-discussed non-twisting case. 



6 Conclusions 

We have presented a systematic study of geometric properties of KS met- 
rics in higher dimensions. We have preliminary discussed general results 
that do not require any further assumptions. Namely, the KS vector k is 
geodetic if and only if the energy-momentum tensor satisfies T a i,k a k b = 
(proposition 1). When this happens (e.g., in vacuum) the Weyl tensor is 
of type II (or more special) and k is a multiple WAND (proposition 2). 
Furthermore, we have shown that optical properties of k are the same with 
respect to both the flat background and the full metric. Subsequently, our 
analysis has focused on vacuum solutions. 

For non-expanding metrics the most general KS solution is now known, 
since they are equivalent to the type N Kundt class in vacuum (propo- 
sition 4). Expanding solutions required a more detailed analysis. Again, 
they turned out to be algebraically special, but the only possible types are 
II and D (proposition 5). We have also shown that the choice of a possi- 
ble KS congruence is restricted in vacuum by an "optical constraint" . In 
n = 4 dimensions this requires k to be a shearfree congruence, in agree- 
ment with the standard Goldberg-Sachs theorem. For n > 4 we have 

non-constant curvature horizons [32], and 'exceptional' solutions with 'zero mass' fi = [32,33] 
(see also [16]). These cannot be KS spacetimes because their Weyl components C'ijki contain 
an r -2 term which is absent from (78). 
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proven a partial and apparently weaker extension of this theorem, which 
however naturally reduces to the familiar result when n = 4. This ex- 
tension has been derived here only for KS solutions but, interestingly, it 
agrees with previous results for general vacuum type III/N spacetimes [8]. 
Moreover, by integration of the Ricci identities we have fixed the depen- 
dence of the optical matrix Lij and of the KS function TC on the affine 
parameter r along k. This enabled us to prove the presence of a curvature 
singularity in "generic" expanding KS spacetimes. As a side remark, let 
us note that we were interested in the n > 4 case, but our analysis also 
applies to the n = 4 case, for which we have rederived various known 
results previously scattered in several publications (many reviewed, how- 
ever, in [4]). 

In future work, it will be interesting to employ our results to possibly 
find new expanding KS solutions, and to understand to what extent the 
Myers-Perry metrics exhaust the expanding vacuum KS class. In addition, 
further investigation will admit non-zero matter fields and generalized KS 
solutions, in which the background is not necessarily flat. These include 
other important spacetimes such as rotating black holes in (A)dS [34, 35] 
and their NUT extensions [25] . 
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A Frame components of the Riemann ten- 



When k is geodesic and afhnely parametrized we find the following frame 
components of the Riemann tensor corresponding to the line element (2) 



iioiii = -6i(DH) + 2L m DH + LjiSjH + 2H{2L Jt L {lj] + L 3l A 3l ), (A5) 



sor 




fioioi = D 2 Ti, Rouj = 2AjiDTi + AUSkyA^, 

Ronj = —LijDTl — 27iAkiLkj , 

Riju = 4Ji.(AijAki + AkyA^i + Si^S^k), 



Roioj — 0, -Roio; — 0, Roijk — 0, 



(Al) 
(A2) 
(A3) 
(A4) 
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+ L lk M(ij) - ZHA ik A jk - L k(i Mj)x - L (i \ k M\j)i). (A7) 

For certain calculations, it may be useful to note that, using the Ricci 
identities (Ilk, [9]), the above component Rujk can also be transformed 
to the somewhat different form 

Ruj k = 2Aj k SiH. — LkiSjU. + LjiSkH + 2H^2S[ k Sj]i + 2 M\j k \ Su 

- 2 Miy S k]l + 2L u A jk - 2L 1[j A k] ^ . (A8) 

Let us finally emphasise that throughout sections 4 and 5 and in ap- 
pendix D we restrict to vacuum spacetimes, so that C a bcd = Rated and 
the Weyl tensor is there given simply by the above expressions. 

B Expanding solutions with DH = do 
not exist in vacuum 

In the special case when DH = 0, from the vacuum equation (33) we get 
uj — 0, i.e. Lij — Sij. Then eq. (30) reads 

Si k Sj k = (n — 2)9Sij. (Bl) 

Using an eigenframe of Sij , it is easy to see that the only possible solution 
has the form 

& i =a (a) diag(l,0,0,...), (B2) 

with S( 2 ) = (n — 2)6. Now, putting DH — 0, Aij — 0, Lij — Sij and 
eq. (B2) into eqs. (A2)-(A4), we find that all components of the Rie- 
mann (Weyl) tensor with non-negative boost weight vanish. The algebraic 
type must thus be III or N. However, by analyzing the Bianchi identities 
it has been proven in [8] that the canonical form of the expansion ma- 
trix Sij for type N and (non-twisting) type III spacetimes in vacuum is 
Sij = s (2) diag(l, 1, 0, . . .) (cf eqs. (50) and (C.20) of [8]). This is clearly 
incompatible with (B2). Therefore, there do not exist vacuum solutions 
of the KS class with 0^=0 and DH = 0. 



C Solving the optical constraint 

In this appendix we provide a solution of the optical constraint (45) for 
the matrix Lij, in the case of an expanding KS vector, i.e. 9 7^ 0. We 
introduce the compact notation 

L lk L lk a 2 + lj 2 + (n - 2)6 2 
(n-2)6 (n-2)6 ' [ } 

so that (45) simply reads 

LuLji = J-Sij. (C2) 

In the following, it will be convenient to analyze separately the two 
possible cases det L / and det L = 0. 
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C.l Non-degenerate case 



Let us start with the non-degenerate case det L 7^ 0, so that there exists 
a matrix inverse of Lij. Let us denote by such inverse matrix (or. 
sometimes, its (y)-element - there will be no ambiguity according to the 
context). It is also convenient to define symbols for its antisymmetric and 
symmetric parts, i.e. 

B H = C'ij = L^y (C3) 

Now, multiplying (C2) by L^- one gets 

L ik =FS ij L k T J 1 , (C4) 
and by further multiplication by L^ 1 we find 8 ik — TCik, so that 

=F- 1 5 ij +B ij . (C5) 

Eq. (C4) thus becomes 

Lik = Sik + TSijBkj- (C6) 

The symmetric and anti-symmetric parts of this equation give rise, 
respectively, to 

BkjSji + BijSjk = 0, (C7) 
2Aik = J- (BkjSji — BijSjk)- (C8) 

The symmetric matrix Sij defines a natural frame of orthonormal 
eigenvectors. It is thus convenient to identify our basis vectors m' 1 ' with 
such eigenvectors, 14 so that Sij = diag(s( 2 ), S( 3 ), S( 4 ), . . .). Eqs. (C7) and 
(C8) then read 

Sfci(a(i)-*(fc)) = 0, (C9) 
2A ik = TB M {s {i) + s (fc) ). (CIO) 

If Bki = for all values of the indices, then eqs. (C9) and (CIO) are 
identically satisfied and, by (C5), Lij — !F5ij, and all eigenvalues of Sij 
take the same value S(j) = T. 

If, instead, for some components B^i / 0, eq. (C9) implies that the 
corresponding eigenvalues of Sij coincide (i.e., s^j = s^))- After ordering 
the basis vectors so as to have multiple eigenvalues next to each other 
on the diagonal of Sij, the matrix Bki must thus be composed of (anti- 
symmetric) blocks with their diagonal on the main diagonal, each block 
corresponding to a repeated S/j^ (e.g., B23 7^ implies S( 3 ) = S( 2 ) and 
then, if S( 4 \ 7^ S( 2 ), we have B24, = 0, etc.). Moreover, there can be at 
most one eigenvalue of Sij with multiplicity 1 , and this must equal T (sup- 
posing, e.g., that S( 2 ) 7^ s^j for k = 3, . . . , n — 1, then one gets B 2 f, = 0, 

14 A priori, there is no reason to expect that this frame is parallely propagated along k, a 
condition that we do not need here. However, we shall see in the following appendix that such 
an eigenframe is, in fact, compatible with parallel transport. 
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so that, by (C5), L^- 1 = = LT 1 , L22 = ^ r_1 , and thus S( 2 ) = T ; and of 
course there cannot be distinct eigenvalues both equal to JF). 

Indices i, j, k, ... take values within a given block (and thus S/j-, = 
S(fc)), so that eq. (CIO) now splits into separate equations for each block 
of the matrix Bki, i.e. 

A rk = s G) FB m , (Cll) 

and components of Aij with indices referring to different blocks necessarily 
vanish (since thus does Bij). The matrix Aij has thus the same block- 
structure of Bij (or simpler, if some of the su\ vanish). Eq. (Cll) also 
implies that all eigenvalues S( 4 ) (for all blocks) must be non-zero, since we 
are now considering the case det L 0. 

Using (C5), (C6), and (Cll), within each block we can thus write 

=F-%- j+ Bjj, Ljj = s G) (5ij + FBjj). (C12) 

But since these two blocks must be inverse to each other (and since 7^ 
takes the same value for any i within a block), we necessarily have 

BrkBjk = T2 W <%■ (C13) 
■r 

This implies that any block-matrix Bij must be even-dimensional (and 
that (T — S(jj)s(j) > 0), unless it vanishes. When this happens, one has 
simply Ljj 1 = J r ~ 1 Sjj and Ljj = s^Sjj, so that 

8&=T (&B rj =0). (C14) 

Hence, if there exist more than one block satisfying Bj-j = 0, to such values 
of indices there corresponds only one possible eigenvalue sm = T of the 
matrix Sij, i.e. Sij contains a diagonal block given by jTdiag(l, 1, . . . , 1). 

Within each even-dimensional block with Bjj 0, = s^Sjj and 
we can thus finally perform a rotation so as to put this anti-symmetric 
matrix in a canonical form with two-dimensional anti-symmetric blocks 
(constrained by (C13)) along the diagonal and zeros elsewhere. 

To summarize, we have shown that the matrix i^ 1 can always be 
written in a block-diagonal form. If its antisymmetric part Bij does not 
vanish identically, there is a number p < (n — 2)/2 of two-dimensional 
blocks of the form 

^U'firr - J' (CI5) 

n - 2 

H=l,...,p, 0<p< 

s (2n) need not be all distinct and thus some two-blocks may be identical. 

In addition, there is a (n — 2 — 2p) x (n — 2 — 2p)-dimensional diagonal 
block 

C' 1 = ^~ 1 diag(l, 1, . . . , 1). (C16) 
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Since each block can be inverted separately, finding the explicit form 
of the matrix Lij is now straightforward. This consists of p blocks of the 
form 



Am) = 



(C17) 



/ r ]i/2' 

S(2 M ) [S(2f,)(^- S(2 M ))J 

-[«(2 M )(-? r - «(2 M ))J «(2 M ) 

and of one (n — 2 — 2p) x (n — 2 — 2p)-dimensional diagonal block 

£ = .Fdiag(l, !,...,!). (C18) 



V 



(n — 2 — 2p) 



More explicitly, the matrix L^ thus takes the form 



C 



(2) 




(C19) 



V 



C.2 Degenerate case 

Let us now study the degenerate case detL = = detL T . This im- 
plies that there exists a non-zero vector v (living in the 'transverse' space 
spanned by vectors m' 1 ') such that L T v = 0. By (C2) this gives also 



Lv = 0. If we now choose an orthonormal basis such that, say, m 



(n-l) 



corresponds to v, we can rewrite L 



J (n-l),i — 



= Lv as 



J i,(n-1) • 



(C20) 



In addition, the optical constraint (C2) now reduces to L-^L-jj. = FSjj, 
with the barred indices ranging from 2 to (n — 2) only. In other words, 
we have the same equation in a reduced space. We can then distinguish 
the two cases det Lj-j and det L^j = 0, and similarly proceed until we 
find a subspace with a non-degenerate reduced matrix L. Finally, we can 
proceed exactly as in the non-degenerate case. It is thus obvious that the 
general form of Lij will have the same block diagonal form (C19), except 
that now the diagonal block contains also (n — 2 — m) zeros, 



C = J r diag(l, ... ,1, 0,... ,0), 

(m — 2p) (n — 2 — m) 



(C21) 



where m = rank(L), and < 2p < m. The non-degenerate case clearly 
corresponds to m = n — 2. 
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D Integration of the Ricci identities 



The general form of the Ricci identities in higher dimensions has been 
given in [9]. Here we consider them only in the case 9 ^ 0, and in a form 
simplified by the fact that I = fe is geodetic and afhnely parametrized. 
Further simplification could be achieved by using a null frame that is 
parallely transported along k (which is always a possible choice [9]), i.e., 

i 

Mjo = = Nio- For certain purposes, however, it is desirable to retain 
the full freedom of null rotations of n, and for now we thus only impose 

M j0 = 0. (Df) 

We also employ the special form that the Riemann tensor takes in 
case of KS spacetimes, as given in appendix A, and use the compact nota- 
tion (62)-(64) (strictly speaking, this is defined for the Weyl tensor, but in 
vacuum we can equivalently apply it to the Riemann tensor). Moreover, 
we keep into account consequences of the vacuum equations discussed in 
appendix C. 

D.l Sachs equations: form of La- 
under the above assumptions, eq. (llg, [9]) becomes 

DLij = —LikLkj- (D2) 



D.l.l Non-degenerate case 



By methods similar to those of [12], the solution to (D2) can be written 
in terms of its inverse as [36] 



L- 1 = rS t] + (L%\ 



(D3) 



where {L {) ) ij 1 represents at r = 0. For compatibility with the op- 
tical constraint, the 'initial value' (L )" 1 must be of the canonical form 
obtained in appendix C. The solution for any r can then be obtained 
immediately by adding a term rSij to eqs. (C15) and (C16) evaluated at 
r = 0. 

Hence, the inverse matrix i^ 1 consists of p two-blocks and of one 
(n — 2 — 2p) x (n — 2 — 2p)-dimensional block of the form, respectively, 



-(m) 



T- U 

~1P 

(2m) 



•^0 



1/2 



~1F 

<2m) 



1/2 



Fo 1 + ■ 



£- 1 = (^ - 1 +r)diag(l,l,...,l). 
Correspondingly, the matrix Lij (C19) consists of blocks 
1 



(D4) 
(D5) 



-(m) 



1 + 2r-6- ( ( ' 2M) + r 2 S ° 2M) ^o 
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^)(l+^0) [4m)(^- S (2m))] 1/2 
1/2 



,(D6) 



° -diag(l, !,...,!). (D7) 



1 + rF 

By construction, this matrix Lij satisfies the optical constraint at 
r — 0. However, note that the last two equations are of the same form of 
(C17) and (C18), respectively, provided one rewrites them using 



^ l + 2r S 2M) +r2 s 2M) ^ ' ' 1 + rT, ' 



This implies that Lij automatically satisfies the optical constraint for 
any r. In other words, the canonical frame determined by the optical 

i 

constraint is compatible with the condition Mjo = (and, in particular, 
with parallel transport along k if N i0 = is also required). 

For practical purposes it will be convenient to simplify the above equa- 
tions by shifting the affine parameter as r — f — l/To, and simultaneously 

defining the new quantities a° 2)l) = [s° 2/j) (To -«(2 M ))] /( s (2 M )-^b)- Then 
one has (dropping tildes over r) 



, *(2/i) 

T -2 . C^O ^2 _„9. 



r a%„, \ ~ 1 



-diag(l, !,...,!). (D9) 



+ K 2M) ) 2 v-«(2m) »■ y 

This general structure of L»j agrees with the results of [16] (up to 
reordering the basis vectors) for the specific case of Myers-Perry black 
holes. 



D.1.2 Degenerate case 

When det L = 0, it is easy to see that the vector v (that satisfies Lv = 0) 
can be parallely transported along k [36]. We can therefore proceed as 
in the non-degenerate case to solve the Sachs equations in the 'non-zero 
block' of Lij, the remaining part of the matrix keeping its zero form. We 
thus do not repeat the results here, just replace the last (n — 2 — m) entries 
of C by zeros in (D9). 



D.1.3 Summary 

In the following, we will often need to refer to the elements of the matrices 
Sij and Aij, and we will need a notation that enables us to handle both 
the non-degenerate and the degenerate case simultaneously. It will also 
be convenient to distinguish between indices a, /3, . . . referring to the non- 
zero eigenvalues of C and p, a, ... to its zero eigenvalues. For further 
compactness, we also use a complex notation (with i being the imaginary 
unit) and we summarize the non-zero elements of Lij as 

r + ia%n) 

S(2 M ) + ^2 m ,2 m +1 = , -, o , (/i = l,...,p) (D10) 
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*(«») = £ (a = 2p + 2,...,m + l), (Dll) 

s (p )=0 (p = m + 2,...,n-l). (D12) 

D.2 Remaining Ricci rotation coefficients 

The above results for Lij enable us to find the r dependence of all Ricci 
rotation coefficients by integrating the corresponding Ricci identities [9] 
that contain a derivative along k. From now on, however, in addition 
to (Dl) we also use the null rotation freedom on n to set 

N i0 = 0, (D13) 

i.e. we use a null frame that is parallely transported along k. 

D.2.1 Form of Lu and Lu 

Let us start from eq. (lib, [9]), which now takes the form 

DLu = -LijLji. (D14) 

These equations decouple according to the block structure of Lij , and the 
resulting coefficients can be compactly written in complex form as 

1-1,2/1 + «il,2 M +l = (Z?,2/i + *'l,2 M +l)( s (2/i) + «-42 M ,2n+l), (D15) 

where Z? j2 /i and Z?,2/i+i are real integration 'constants', and S( 2m ) and 
-42 M ,2 M +i are given in (D10). 

The solution for coefficients Lu with i corresponding to the diagonal 
block of Lij is simply 

Lia = — , L lp = l° lp . (D16) 
r 

Ricci equations (lie, [9]), i.e. 

DLn — —LtjLji, (D17) 
lead to similar expressions for the coefficients Lu, namely 

^2/1,1 + 1^2/1+1,1 = (^2/1,1 + ^2/1+1,1 )( s (2/i) — ^2/1,2/1+1 ), (D18) 

I a i = ^, L pl = l° pl . (D19) 
D.2. 2 Form of M jk 

i 

Next, the coefficients Mjk are analogously determined by eq. (lln, [9]), 
which here takes the form 

DM jk = -MjiLi k . (D20) 

Again, the solutions can be given in natural pairs, each for each block of 
Lij, as distinguished by the value of the last index k, i.e. 

i i j (J i 

Mj, 2 M +lMj, 2 M+ 1 = ("lj,2/i + i m j,2 M +l )(«(2/i) + iMp.an+i), (D21) 
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with real integration 'constants'. When k corresponds to the diagonal 
part of Lij we have simply 



i 
TTl j a 



Mj P = rrijp. 



(D22) 



i i q j o 

Because of the index symmetries of Mjk [8], we require m jfe + m ife =0 
for any i,j,k = 2,...,n — l. 



D.2.3 Form of Mji 

Eq. (11m, [9]) here becomes 

DMj! = -M jk L kl - 2$4. (D23) 
From (D10), one has = -2s (i )Aij and thus, by (A2), $^ = -D(HAij). 

i 

Using also previous results for Mjk we find 



Mji= J2 



+ 



r i o iO . 

( m i,2M '2/i,l + m J,2M+l '2 M +1,1 J 



- 2 + K 2 ,)) 2 



(2m) c * i() * i \ 

—1 ^2(, m j,2 M «2 M +l,l - m 3,2n+l ( 2 M ,lJ 

(2/0 



i o 

+ — — — - m JP + 2A l3 U + mji , 



(D24) 



j Q j 

where the real 'constants' satisfy rriji + mn = 0. 



D.2.4 Form of L n 

One can also easily integrate (11a, [9]), here reduced to 

DLn — —L\iLi\ — i?oioi, 
to fix the r-dependence of Ln, i.e. 

p r 

ill = — -DH + l"x + I ('l, 2^2^,1 + ^l,2 M +l'2 M +l,l)S(2 M ) 



|tx=l 



('l, 2^^2^+1,1 — 'l,2fj+1^2M,l)^42M,2M+lj 



;0 ,0 

L la l al ,0 ,0 
l lp L pl r - 



(D25) 



(D26) 



D.2.5 Form of N, 



Eq. (llj, [9]) here reads 



DN jk = -NjiUk - $kj. 



(D27) 
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Note that, for k = 2/i, 2/i + 1, it is more natural to deal with a complex 
unknown, namely Nj^y. + iNj^+i ■ Then, using (62), (63) and (D10)- 
(D12), eq. (D27) can be solved obtaining 

Nj,2f, + iNj t 2p+l = ("-3,2m + * n j,2M+l) ( S (2m) + ^2 M ,2 M +l) 

+ W(S(2 M ) - M 2M ,2 M +l)((5j,2 M + i6j,2n+l), (D28) 

N ja = \ (n% + HS ja ) , N jp = n° jp , (D29) 
where n^- are real 'constants'. 

D.2.6 Form of Nn in an adapted parallely transported 
frame 

So far, the frame vector n has not been specified, except for the require- 
ment that it be parallely transported along fe. While retaining the latter 
condition, we can still perform a null rotation 

fe = fe, n = n + Zi m {i) - ^z k z k £, m (i) = m (i) - zd, (D30) 

provided Dzi — 0, so as to simplify some of the Ricci rotation coefficients. 
In particular, if we choose 

Zi = -la (for i = 2, . . . , m + 1), z p = 0, (D31) 

we obtain (see the transformation properties given in [9] - we drop hats 
from the transformed coefficients) 

L a =0 (i = 2,...,m+l), (D32) 

which is equivalent to Z^.i = ^+1,1 = C = 0, while the L p \ — l pl are 
unchanged. (Alternatively, one can also set Lu = 0, but this appears to 
be less convenient in what follows.) 

With this choice, eq. (1 If, [9]) becomes 

DN n = -n° ip l° pl + R 101i . (D33) 

In addition, using eq. (D14) and the commutator (22, [21]), eq. (A5) 
reduces to 

RoiH = -D(5iH)-2D(L u H) (i = 2, . . . , m + 1), (D34) 
Roup = ~D(S P H) - 2D(L lp H) + l° pl DH. (D35) 

We can thus straightforwardly integrate eq. (D33) and find 

Nn = -n°Xir + Sm + 2LuH + n° n , (i = 2,...,m + 1) (D36) 
N P i = -n° p Xir + 8 P H + (21% - l° pl )H + n° pl , (D37) 

where Tt and Lu are specified by (56) and (D15), (D16). For certain 
applications it may be useful to recall that the coefficients Nn vanish if and 
only if the frame vector n is geodetic (an affine parameter corresponding 
to in = 0). 
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Note finally that, when l° Pi i = '2^+1,1 = Si = 0, eq. (D24) simplifies 

to 

Mji = - rn° p l° pl r + 2A i3 U + m°i , (D38) 
and eq. (D26) becomes 

L11 = -DH + 2?i - l° p l° pl r. (D39) 
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